Adjustment of control parameters of section of optical fiber network

ABSTRACT

Adjustment of one or more control parameters of a section of an optical fiber network involves taking measurements of optical signals in the section, deriving estimated data from the measurements and from knowledge of the section, where the estimated data is a function of optical nonlinearity and of amplified spontaneous emission, and applying one or more control algorithms using the estimated data to adjust the one or more control parameters.

TECHNICAL FIELD

This document relates to the technical field of optical communications and specifically to the control of components in an optical fiber network.

BACKGROUND

Current best practices for determining optical parameters in an optical fiber network look at equalizing the ratio of amplified spontaneous emission (ASE) to signal power on channels over an optical section while respecting channel power limits to manage the fiber optical nonlinear effects. This equalization addresses the strong power tilt that can accumulate across spans of optical fiber mainly due to Stimulated Raman Scattering (SRS). These methods rely heavily on offline simulations to determine good control parameters, such as peak power. This is operationally burdensome and error prone.

U.S. Pat. No. 9,438,369 describes increasing capacity by optimization after nonlinear modeling. U.S. Pat. No. 8,364,036 describes controlling optical power within domains, and exchanging state information between domains. U.S. Pat. No. 8,781,317 describes methods to measure phase nonlinearities. U.S. Pat. No. 7,894,721 describes global optical control where receiver changes are correlated to network perturbations. U.S. Pat. No. 7,457,538 describes performance monitoring using the analog-to-digital converter of the receiver. U.S. Pat. No. 7,376,358 describes location-specific monitoring of nonlinearities. U.S. Pat. No. 7,356,256 describes digital monitoring along the optical line. US Patent Publication No. 2016/0315711 describes controlling the optical spectral density in a section.

SUMMARY

Through the latest innovations, optical networks are capable of dynamically changing optical paths, and flexible transceivers are capable of changing modulation formats and other transmission parameters. In this environment, optical line control that provides good performance, scalability, and self-optimization is desirable.

Adjustment of one or more control parameters of a section of an optical fiber network involves taking measurements of optical signals in the section, deriving estimated data from the measurements and from knowledge of the section, where the estimated data is a function of optical nonlinearity and of amplified spontaneous emission, and applying one or more control algorithms using the estimated data to adjust the one or more control parameters.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 illustrates a method for adjustment of control parameters in section of an optical fiber network;

FIG. 2 illustrates an example section of an optical fiber network;

FIG. 3 illustrates an example concave value function of excess margin; and

FIG. 4 illustrates a first derivative of the example concave value function.

DETAILED DESCRIPTION

Optical network topologies can range from simple unamplified point-to-point, to branching chains of reconfigurable optical add drop multiplexer (ROADM) sections, up to a full multi-connected mesh that spans a continent.

In wavelength division multiplexing (WDM) systems, an optical fiber network connects wavelength selective switch (WSS) nodes via spans of optical fibers and optical amplifier devices. Pairs of flexible coherent transceivers are connected over paths through the optical fiber network. Different channels are propagated through different paths in the network. A flexible coherent transceiver can be reconfigured allowing transmission parameters, for example, modulation scheme, to be adapted to the selected path.

Some elements of the optical fiber network have the ability to do some level of per-channel power control. Such elements may include, for example, the transmitter portions of the flexible coherent transceivers, a variable optical attenuator (VOA) under software control, and optical equalizers. In another example, per-channel power is controllable by provisioning a wavelength selective switch (WSS) node with loss values. A WSS node has switching capabilities and per-channel power control.

Some elements of the optical fiber network have the ability to do some level of total power control. Such elements may include, for example, optical amplifier devices. For example, the gain of an optical amplifier device is controllable by provisioning the optical amplifier device with a target gain. Equivalently, the total output power (TOP) of an optical amplifier device is controllable by provisioning the optical amplifier device with a target total output power.

Some optical amplifier devices also have the ability to do some level of per-channel power control, by provisioning the optical amplifier device with a target gain tilt profile. For simplicity, this document focuses on the following control parameters of a section of an optical fiber network: the loss values of a WSS node, which affect the launch powers of the signals launched into the optical fibers, and the target gain values (or target TOP values) of optical amplifier devices.

FIG. 1 is a flowchart illustration of a method for adjustment of control parameters in a section of an optical fiber network. A section may comprise most or all of the optical fiber network. If the optical fiber network is small, the section may indeed comprise all of the network. However, it is generally advantageous for the method to control a single point-to-point section of optical amplifier devices and spans of optical fiber between two nodes that contain ROADM, WSS, or other switching hardware that may be present.

At 2, measurements of optical signals are taken at various locations in the section. The measurements may include per-channel optical power (also referred to as power spectral density, especially in a flexible grid system) and total output power. For example, an optical power monitor (OPM) device is able to measure per-channel optical power by switching the optical connection to its input. Due to the cost of an OPM device, there is generally not an OPM device at each optical amplifier device. Taps and photodiodes may be placed, for example, at the input and at the output of the optical amplifier devices. Each photodiode is operative to measure the total optical power at the location of the tap. At locations where there is an OPM device and a tap and photodiode, the measurement of total optical power may be used to calibrate the per-channel optical power measured by the OPM device.

At 4, estimated data is derived from the measurements, from the target values, and from knowledge of the section and its components. The estimated data may include, for example, the estimated per-channel optical power at the output of the optical amplifier devices, the estimated incremental amplified spontaneous emission (ASE) power introduced by the optical amplifier devices, and the estimated self-phase modulation (SPM) and cross-phase modulation (XPM) variance accumulated in the section. The estimated data may be derived using a modeling engine that models the propagation of signals through the components of the section. Alternatively, the estimation of nonlinearities and noise may be derived from specific measurements of parameters as described in U.S. Pat. No. 8,594,499, U.S. Pat. No. 7,356,256, U.S. Pat. No. 6,128,111, U.S. Pat. No. 6,687,464, U.S. Pat. No. 6,839,523, U.S. Pat. No. 7,376,358, U.S. Pat. No. 6,072,614, U.S. Pat. No. 6,064,501, and U.S. Pat. No. 5,513,029.

The estimated data is then used in a control algorithm to adjust the control parameters. Various control algorithms are contemplated. For example, the control algorithm may make use of gradients and slew-rate limited steepest descent. At 6, gradients of an objective function are evaluated, using the measurements and the estimated data. The values of the gradients are inaccurate, for at least the reason that the measurements are noisy, the knowledge of the section and its components may be inaccurate or incomplete, the modeling engine is inaccurate, and the estimated data is inaccurate. Some of the channels propagated through the section carry live traffic. That is, some of the channels are in-service channels carrying traffic for customers. It is therefore important not to adjust the components of the section in a manner that would jeopardize or degrade or destabilize the in-service channels.

At 8, the values of the gradients are used in steepest descent algorithms to adjust control parameters of the section by a small step in a direction of optimization of the objective function. That is, small adjustments are applied to target values such as loss values of a WSS node and the target gain (or target total output power) of an optical amplifier device. Steepest descent algorithms are known to be noise tolerant and to be very safe if small steps are taken. The values of some control parameters that are adjusted may be set points for algorithms that control other control parameters. For example, a value of a per-channel optical power out of a WSS node may be a set point for an algorithm that adjusts the loss of the relevant pixels of that WSS node. A total power may be a set point for an algorithm that adjusts total gain, which may be a set point for a digital control loop which adjusts a VOA loss and pump currents. A pump current may be a set point for an analog loop which adjusts a Field Effect Transistor (FET) bias.

The method illustrated in FIG. 1 may be repeated over the lifetime of use of the optical fiber network. For example, the method may be repeated every few seconds for 25 years. It is not necessary that all control parameters be adjusted in each iteration of the method. Various changes occur over time, yielding updated measurement data, updated estimated data, updated values for the gradients, an updated direction of optimization of the objective function, and updated adjustments to the control parameters.

The optical fiber network may be partitioned into sections arbitrarily. For simplicity, this document focuses on an example section that enables transmission of a set of optical signals along a particular transmission direction from a first WSS node to a second WSS node. (Signals are also directed along the opposite transmission direction, where the roles of ingress and egress are reversed. However, so as not to obscure the description of the technology, transmission along that opposite direction is not illustrated and is not discussed in this document.)

FIG. 2 illustrates an example section 10 of an optical fiber network. An ingress WSS node 12 is connected to an egress WSS node 14 via spans 16 of optical fiber. The length of a span 16 of optical fiber is typically in the range of approximately 80 km to approximately 100 km. The spans 16 of optical fiber are coupled via optical amplifier devices 18. An optical pre-amplifier device 20 in the ingress WSS node 12 is optically coupled to the first span 16 of optical fiber. An optical pre-amplifier device 20 in the egress WSS node 14 is optically coupled to the final span 16 of optical fiber. One can index the spans 16 and the optical (pre-)amplifier devices 18,20 by an index j, with N representing the total number of spans of optical fiber coupling the ingress WSS node 12 to the egress WSS node 14. For example, one can refer to the optical pre-amplifier device 20 in the ingress WSS node 12 as the first optical amplifier device, whose output is launched into the first span of optical fiber. Similarly, the output of the optical amplifier device j is launched into the span j of optical fiber.

As discussed above, the measurements taken at various locations in the section may include per-channel optical power (also referred to as power spectral density, especially in a flexible grid system) measured by OPM devices and total output power measured by photodiodes. In the example section 10, OPM devices 22 at the ingress WSS node 12 and at the egress WSS node 14 are able to measure per-channel optical power across the spectrum at the output of the respective optical pre-amplifier device 20. In the example section 10, taps and photodiodes are present at the input and at the output of each optical (pre-)amplifier device 18,20 and are illustrated in FIG. 2 by small black squares. Each optical amplifier device 18 is comprised, together with its respective taps and photodiodes and together with a shelf processor 24, in a network element 26. For simplicity, only one such network element 26 is illustrated in FIG. 2. There is a shelf processor 28 comprised in the ingress WSS node 12 and a shelf processor 30 comprised in the egress WSS node 14.

Each photodiode is operative to measure the total optical power at the location of its respective tap. At the output of the optical amplifier device j, the photodiode measures the total output power, which includes both optical signal power and ASE power. The per-channel power measured by the OPM device 22 is reliable only in terms of relative power across the spectrum, because the loss along a cable 32 coupling the output of the optical pre-amplifier device 20 to the OPM device 22 is generally not known. The total output power measured by the photodiode at the output of the optical pre-amplifier device 20 in the ingress WSS node 12 can be used to calibrate the per-channel optical power measured by the OPM device 22, thus yielding a calibrated set of per-channel optical power measurements {P₁[i]}, where P₁[i] is the power of the channel i launched into the first span of optical fiber. The total output power measured by the photodiode at the output of the optical pre-amplifier device 20 in the egress WSS node 14 can be used to calibrate the per-channel optical power measured by the OPM device 22, thus yielding a calibrated set of per-channel optical power measurements {P_(N+1)[i]}, where P_(N+1)[i] is the power of the channel i output from the optical pre-amplifier device 20 in the egress WSS node 14. The integration of per-channel optical power measurements to yield an aggregate power, comparison of the aggregate power to the measured total optical power, and calibration may be performed by firmware (not shown) in the WSS node 12,14. Alternatively, the integration, comparison and calibration may be performed by any other suitable firmware executed by a processor within the example section 10. Conventional optical power units are dBm. In this document, the per-channel optical power measurements {P_(j)[i]} are conveniently measured in units of Nepers relative to a Watt, because it is more convenient for the calculus of equations appearing hereinbelow.

A control system embedded in the section is operative to provision certain components of the section with specific target values. For example, the control system is operative to provision the ingress WSS node 12 with loss values, and to provision the optical amplifier devices 18 and the optical pre-amplifier devices 20 with respective target gain values or target TOP values. The control system comprises, for example, hardware (not shown) located in the ingress WSS node 12, hardware (not shown) located in the optical amplifier devices 18 and in the optical pre-amplifier devices 20, and control firmware 34 executed by any one of the shelf processors within the section 10, for example, the shelf processor 30 comprised in the egress WSS node 14. The control firmware 34 is stored in non-transitory computer-readable media coupled to the shelf processor.

In an alternative implementation, the control firmware 34 is executed by an external processor (not shown) that is in communication with the controllable elements of the section. The external processor may be located in a physical server or may be virtualized as part of a cloud infrastructure. The apparatus in which the external processor is located may also store the control firmware 34 in non-transitory computer-readable media that is accessible by the external processor.

As discussed above, estimated data is derived from the measurements, from the target values, and from knowledge of the section and its components. The estimated data may be derived using a modeling engine that models the propagation of signals through the components of the section.

The knowledge of the section and its components may include “known characteristics”. Manufacturers and/or distributors of the components may provide some of the known characteristics. Other known characteristics may be determined by testing and/or calibrating the components. Still other known characteristics may be provided by inspection of the section. The known characteristics may include, for example, the topology of the section, one or more optical amplifier characteristics such as amplifier type (e.g. Erbium-doped fiber amplifier (EDFA), distributed Raman amplifier, lumped Raman amplifier), noise figure, ripple, spectral hole burning, and Total Output Power (TOP) limit, and one or more optical fiber characteristics such as fiber type, span length, nonlinear coefficients, effective area, loss coefficients, total loss, chromatic dispersion, and Stimulated Raman Scattering (SRS).

Measured data (raw and/or calibrated), control data, and (optionally) known characteristics, are communicated within the section 10 over an optical service channel (OSC), also known as an optical supervisory channel. The WSS nodes 12,14 and the network elements 26 each comprise circuitry 36 to support the OSC.

The modeling engine models the propagation of signals through components of the section. Specifically, the modeling engine employs fiber models for the spans 16 of optical fiber in the section 10 and employs amplifier models for the optical (pre-)amplifier devices 18,20. Modeling firmware 38 that uses the modeling engine is executed by any one of the shelf processors within the section 10, for example, the shelf processor 24 comprised in the network element 26. The estimated data derived by the modeling engine may include, for example, the estimated per-channel optical power {P_(j)[i]} at the output of the optical amplifier j, where P_(j)[i] is the power, measured in units of Nepers, of the channel i launched into the span j of optical fiber, and the estimated incremental ASE power {ASE_(j)[i]} at the output of the optical amplifier j. The modeling engine may employ known techniques to derive the power evolution of the optical signals through the section and to derive the incremental ASE power.

The accuracy of the estimated per-channel optical power at each of the fiber interfaces is important. Stimulated Raman Scattering (SRS) may impart in the range of approximately 1 dB to approximately 2 dB power tilt across the C band (1525 nm to 1565 nm) and in the range of approximately 3 dB to approximately 4 dB power tilt across the L band (1565 nm to 1610 nm). These power tilts may accumulate between spans where there is no WSS node to equalize the tilts. Channels at different optical powers experience very different optical degradation in terms of ASE (at low channel power) and optical nonlinearities (at high channel power). Good modeling of the SRS tilt per span of optical fiber is part of what contributes to accurate estimated per-channel optical powers and accurate estimated incremental ASE powers.

Once the estimated per-channel optical powers are of sufficient accuracy (which could be determined, for example, by comparing the estimated per-channel optical powers for the output of the optical pre-amplifier device 20 in the egress WSS node 14 with the calibrated set of per-channel optical power measurements {P_(N+1) [i]}), the modeling engine may derive the estimated self-phase modulation (SPM) and cross-phase modulation (XPM) variance accumulated in the section 10. The estimated data is thus a function of optical nonlinearity and of ASE.

The modeling engine may model nonlinear interactions within the spans 16 of optical fiber as Gaussian noise, as described in P. Poggiolini, “The GN Model of Non-Linear Propagation in Uncompensated Coherent Optical Systems”, Journal of Lightwave Technology, Vol. 30, No. 24, Dec. 15, 2012; P. Poggiolini et al. “The GN Model of Fiber Non-Linear Propagation and its Applications”, Journal of Lightwave Technology, Vol. 32, No. 4, Feb. 14, 2014. Alternatively, the modeling engine may employ a different model of the nonlinear interactions, for example, full non-linear Schrodinger Equation solutions using Fast Fourier transform (FFT) or finite difference methods.

As described above, gradients of an objective function are evaluated, using the measurements and the estimated data.

In one aspect, the goal of the objective function is to minimize the total degradation through the section. Optimization of this objective function minimizes a weighted sum of ratios of the total noise power from ASE and optical nonlinearities to the power of the optical signals. This objective function is suitable for systems where there is no software connection to convey information from the receiver modem to the section.

An example objective function V₁ for a section, with the goal of minimizing the total degradation through the section, is given in Equations (1) and (2):

$\begin{matrix} {V_{1} = {{\sum\limits_{i = 1}^{N_{CH}}{{C\lbrack i\rbrack}{{LNSR}\lbrack i\rbrack}}} = {\sum\limits_{i = 1}^{N_{CH}}{\sum\limits_{j = 1}^{N}{{C\lbrack i\rbrack}{{LNSR}_{j}\lbrack i\rbrack}}}}}} & (1) \\ {{{LNSR}_{j}\lbrack i\rbrack} = {\frac{A\; S\;{E_{j}\lbrack i\rbrack}}{e^{P_{j}{\lbrack i\rbrack}}} + {\sum\limits_{k = 1}^{N_{CH}}{{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2{P_{j}{\lbrack k\rbrack}}}}}}} & (2) \end{matrix}$

In Equation (1), LNSR [i] denotes the incremental line noise-to-signal ratio (NSR) for the channel i in the section, which can be expressed as a weighted sum over spans j of optical fiber of the incremental line NSR for the channel i in the span j, denoted LNSR_(j)[i]. C[i] is a weighting value for the channel i to optionally bias the objective function for particular higher-value signals. N_(CH) denotes the number of channels in the signals in the section, and N denotes the number of spans j of optical fiber in the section. C[i] may be a customer-defined weighting value. Alternatively, C[i] may be defined in a different manner. For example, when C[i]=Baudrate [i] the objective function V₁ will converge to control parameters that maximize the capacity-bandwidth product in the optical fiber network. In another example, when C[i]=Baudrate [i]×SNR[i] where SNR[i] is an estimate of the signal-to-noise ratio (SNR) in linear units at the receiver modem whose channel i traverses the section, the objective function V₁ will converge to control parameters that maximize the capacity of the optical fiber network.

In Equation (2), P_(j)[i] is the power of the channel i at the output of the optical amplifier j, which is launched into the optical fiber of the span j, ASE_(j)[i] is the incremental ASE power on the channel i at the output of the optical amplifier j, and NL_(j)[i,k] is the SPM/XPM nonlinear coefficient for Kerr interactions between the channel i and the channel k at the span j. The power P_(j)[i] is measured in units of Nepers relative to a Watt. The second term in Equation (2) is a summation over all channel powers that impact the LNSR of the channel i at the span j. Where nonlinear interactions in one span are independent of nonlinear interactions in another span, the local optimum of this objective function V₁ is the global optimum. The paper I. Roberts, J. M. Kahn, D. Boertjes, “Convex Channel Power Optimization in Nonlinear WDM Systems using Gaussian Noise Model”, Journal of Lightwave Technology, Vol. 34, No. 13, Jul. 1, 2016 proves that the second term in Equation (2) is a convex function in the power P₁ [i] when assuming a Gaussian noise nonlinearity model.

The following discussion derives the gradients of the example objective function V₁, which are evaluated to provide a direction for adjustment of control parameters. A gradient vector ∇V_(1j) for control of the ingress WSS node 12 is derived. A gain gradient for TOP control of the optical amplifier devices 18 is derived.

The example objective function V₁ given in Equation (1) can be expressed as the sum over all spans j in the section of an example span objective function V_(1j), which is given in Equation (3): V _(1j)=Σ_(i=1) ^(N) ^(CH) C[i]LNSR_(j) [i]  (3)

The partial derivative of the example span objective function V_(1j) with respect to channel power for channel i in the span j of optical fiber is given by Equation (4):

$\begin{matrix} {\frac{\partial V_{1j}}{\partial{P_{j}\lbrack i\rbrack}} = {{\sum\limits_{k = 1}^{N_{CH}}{{C\lbrack k\rbrack}\frac{\partial{{LNSR}_{j}\lbrack k\rbrack}}{\partial{P_{j}\lbrack i\rbrack}}}} = {{{- {C\lbrack i\rbrack}}\frac{A\; S\;{E_{j}\lbrack i\rbrack}}{e^{P_{j}{\lbrack i\rbrack}}}} + {\sum\limits_{k = 1}^{N_{CH}}{2{C\lbrack k\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2{P_{j}{\lbrack i\rbrack}}}}}}}} & (4) \end{matrix}$

where the channel power P_(j)[i] is fixed and the sum is over XPM/SPM terms over all channels dependent on channel power P_(j)[i].

A gradient vector ∇V_(1j) for a span j comprises the partial derivative

$\frac{\partial V_{1j}}{\partial{P_{j}\lbrack i\rbrack}}$ for each channel i from 1 to N_(CH). The partial derivative of the example objective function V₁ with respect to WSS loss for the channel k in the section is given by Equation (5):

$\begin{matrix} {\frac{\partial V_{1}}{\partial{P\lbrack k\rbrack}} = {{\sum\limits_{j = 1}^{N}\frac{\partial V_{1j}}{\partial{P_{j}\lbrack k\rbrack}}} = {\sum\limits_{j = 1}^{N}\left\lbrack {{{- {C\lbrack k\rbrack}}\;\frac{{ASE}_{j}\lbrack k\rbrack}{e^{P_{j}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{C\lbrack k\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}}} & (5) \end{matrix}$ where P[k] is the power of the channel k out of the ingress WSS node 12 which affects all spans in the section.

A gradient vector ∇V₁ for the section comprises the partial derivative ∂V₁/∂P[k] for each channel k from 1 to N_(CH). The gradient vector ∇V₁ can be evaluated from the customer values C[i], the known characteristics NL_(j)[i,k], and the measured or estimated data {P_(j)[i]} and {ASE_(j)[i]}.

The incremental ASE power on the channel i induced by the optical amplifier device j is given by the well known Equation (6): ASE_(j) [i]=h*v*B _(e)(NF _(j) [i]*G _(j) [i]−1)  (6) where h is Planck's constant, v is the optical frequency, B_(e) is the electrical bandwidth of the noise filtering in the receiver, NF_(j)[i] is the noise figure for the channel i of the optical amplifier device j, and G_(j)[i] is the gain for the channel i of the optical amplifier device j.

The gradient of the ratio of the ASE power to the signal power term with respect to gain G_(j)[i] is given by Equation (7):

$\begin{matrix} {{\bigtriangledown_{G_{j}{\lbrack i\rbrack}}\left( \frac{{ASE}_{j}\lbrack i\rbrack}{e^{P_{j}{\lbrack i\rbrack}}} \right)} \approx {- \frac{h \cdot v \cdot {{NF}_{j}\lbrack i\rbrack} \cdot B_{e}}{e^{P_{j + 1}^{IN}{\lbrack i\rbrack}}}}} & (7) \end{matrix}$ where P_(j+1) ^(IN)[i] is the channel power at the input to the next optical amplifier device. This gradient is approximately the negative of the ratio of the incremental ASE power to signal power of the next optical amplifier device.

The gain gradient for the optical amplifier device j, averaged over all wavelengths, can be formed as given by Equation (8):

$\begin{matrix} {\frac{\partial V_{1j}}{\partial G_{j}} = {\sum\limits_{k = 1}^{N_{CH}}\left\lbrack {{{- {C\lbrack k\rbrack}}\frac{{ASE}_{j + 1}\lbrack k\rbrack}{e^{P_{j + 1}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{C\lbrack i\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}} & (8) \end{matrix}$

The gain gradient for the optical amplifier device j, averaged over all wavelengths, can be evaluated from the customer values C[i], the known characteristics NL_(j)[i,k], and the measured or estimated data {P_(j)[i]} and {ASE_(j)[i]}.

As described above, the values of the gradients are used in steepest descent algorithms to adjust control parameters of the section by a small step in a direction of optimization of the objective function.

Small adjustments are applied to loss values of a WSS node and to target TOP values of optical amplifier devices. The steepest descent algorithm is applied to the WSS node while assuming that the gains of the optical amplifier devices are fixed. The steepest descent algorithm is applied to all of the optical amplifier devices in parallel while assuming that the WSS pixel drive values are fixed.

For example, two loops may be run in parallel with a decoupling factor, as expressed in the vector Equation (9), Equation (10) and Equation (11):

$\begin{matrix} {{WSS\_ PowerTarget}_{NEW} = {{WSS\_ PowerTarget} - \left( {\frac{MAXSTEP}{\max\left( {\bigtriangledown\; V_{1}} \right)}*\left\lbrack {{\bigtriangledown\; V_{1}} - {{mean}\left( {\bigtriangledown\; V_{1}} \right)}} \right\rbrack} \right)}} & (9) \\ {{TOP\_ Target}_{NEW} = {{TOP\_ Target} - {{{sign}\left\lbrack \frac{\partial V_{1j}}{\partial G_{j}} \right\rbrack}*0.1*{MAXSTEP}}}} & (10) \\ {{{{if}\mspace{14mu}{TOP\_ Target}_{NEW}} \geq {TOP}_{LIMIT}},{{{set}\mspace{14mu}{TOP\_ Target}_{NEW}} = {TOP}_{LIMIT}}} & (11) \end{matrix}$ where WSS_PowerTarget_(NEW) and WSS_PowerTarget have values for each channel k from 1 to N_(CH), the decoupling factor in this example is 0.1, and the target TOP for the optical amplifier device j is subject to an upper limit. An example MAXSTEP is 0.2 dB.

TOP control is used to decouple incremental SNR optimization in this section from changes occurring in other sections of the optical fiber network. The subtraction of the change in average power (which is denoted mean (∇V₁) in Equation (9) but is not quite equal to the average of the changes) keeps the WSS output power constant, and the power launched into the first span is controlled by the TOP of the first amplifier. This scaling is important when one or more of the amplifiers reaches the limit of their TOP and can provide no more power. With this scaling, the allocation of that limited power between the wavelengths is cleanly optimized.

Note also that there is no reliance on a second derivative for step size, and that this simple algorithm is robust to noise.

As mentioned above, where nonlinear interactions in one span are independent of nonlinear interactions in another span, the local optimum of this objective function V₁ is the global optimum. Operationally, this permits the adjustment of the control parameters for one section to be performed in parallel to the adjustment of the control parameters for other sections of the optical fiber network. For example, the two loops expressed in the vector Equation (9), Equation (10) and Equation (11) may be run in parallel independently for several sections of the optical fiber network.

In another aspect, the goal of the objective function is to maximize the capacity or the reliability or both of a network by allocating margin to channels that are at higher risk of failure at their designated capacities by taking away margin from channels with plenty of margin. This objective function is suitable for systems where, for at least some channels, there is a software connection to convey information to the section (or to the external processor) from the receiver modem that receives that channel. This objective function can also be used protect channels in service while trialing a new channel to see if it can sustain a particular high capacity. Another value of this objective function is to assist channels that are experiencing a slow low-probability degradation event such as polarization dependent loss (PDL) by improving this weakened channel's line SNR at the expense of other channels that have higher margin.

An arbitrary concave value function ƒ is introduced that takes as its argument the excess margin SNR_(M)[i] on the channel i as determined at the receiver modem. A positive value for SNR_(M)[i] indicates that total SNR (including ASE, nonlinear effects, and internal receiver modem noise) currently experienced by the channel i exceeds the SNR required for error-free communications on that channel. A negative value for SNR_(M)[i] indicates that the total SNR currently experienced by the channel i is less than the SNR required for error-free communications on that channel. The concave value function ƒ(SNR_(M)[i]) expresses the utility of extra margin on a channel and whether the channel is better off sharing its excess margin. FIG. 3 illustrates an example concave value function ƒ having desirable properties, and FIG. 4 illustrates a first derivative ƒ¹ of the example concave value function. It is scaled so that ƒ(0)=0 and ƒ¹(0)=1. The example concave value function ƒ is neutral (has a value of zero) for zero excess margin, decreases rapidly for negative excess margin, and increases then quickly plateaus for positive excess margin.

An example objective function V₂ for a section that incorporates information from the receiver modem is the sum over all controllable channels of a concave value function ƒ of the excess margin, as given in Equations (12) and (13):

$\begin{matrix} {V_{2} = {\sum\limits_{i = 1}^{N_{CH}}{{C\lbrack i\rbrack}{D\lbrack i\rbrack}{f\left( {{SNR}_{M}\lbrack i\rbrack} \right)}}}} & (12) \\ {{{SNR}_{M}\lbrack i\rbrack} = {{- 10}\;{\log\left( \frac{{LNSR}_{M}\lbrack i\rbrack}{{BLNSR}_{M}\lbrack i\rbrack} \right)}}} & (13) \end{matrix}$

In Equation (12), C[i] is a customer-defined weighting value for the channel i to optionally bias the objective function for particular higher-value signals. D[i] is a metric that is a function of the geographic distance travelled by the channel i from the transmitter to the receiver, or other such network value. Adjusting the function for D allows an adaptation of the trade-off between the use of an optical fiber (if optical fiber on this route is a scarce resource and installing or acquiring rights to more would be very expensive, then, for example, set D[i]=1), and minimizing the cost of the transceivers (if optical fiber is plentiful, then, for example, set D[i]=distance[i]).

In Equation (13), LNSR_(M)[i] is the line NSR for the channel i as measured at the receiver modem, and BLNSR_(M)[i] is a budgeted line NSR which factors in margin, implementation noise, and target Required Noise to Signal Ratio (RNSR), required for the modem to be error free under nominal conditions, from the capacity commitment and forward error channel (FEC) performance for the channel i.

There are many different ways in which the budgeted line NSR for the channel i, BLNSR_(M)[i], can be defined. For example, the budgeted line NSR may be defined as given in Equation (14):

$\begin{matrix} {{{BLNSR}_{M}\lbrack i\rbrack} = {\frac{1}{m_{P}}\left( {{{FEC\_ NSR}\lbrack i\rbrack} - {{INSR}\lbrack i\rbrack} - m_{A}} \right)}} & (14) \end{matrix}$

In Equation (14), FEC_NSR[i] is the NSR for the modulation format for the channel i at the FEC threshold, INSR[i] is the modem implementation noise for the channel i, m_(P) is a multiplicative margin applied to the line NSR, and m_(A) is an additive noise margin.

The example objective function V₂ given in Equation (12) can be expressed as the sum over all channels i of an example channel objective function V₂[i], which is given in Equation (15): V ₂ [i]=C[i]D[i]ƒ(SNR_(M) [i])  (15)

The partial derivative of the example channel objective function V₂[i] with respect to channel power for channel i in the span j of optical fiber is given by Equations (16), (17) and (18):

$\begin{matrix} {\frac{\partial{V_{2}\lbrack i\rbrack}}{\partial{P_{j}\lbrack i\rbrack}} = {{{A\lbrack i\rbrack}\frac{\partial}{\partial{P_{j}\lbrack i\rbrack}}\left( {\sum\limits_{j = 1}^{N}{{LNSR}_{j}\lbrack i\rbrack}} \right)} = {{A\lbrack i\rbrack}\frac{\partial{{LNSR}_{j}\lbrack i\rbrack}}{\partial{P_{j}\lbrack i\rbrack}}}}} & (16) \\ {{A\lbrack i\rbrack} = \frac{{- 10}\;{C\lbrack i\rbrack}{D\lbrack i\rbrack}{f^{1}\left( {{SNR}_{M}\lbrack i\rbrack} \right)}}{2.30\;{{LNSR}_{M}\lbrack i\rbrack}}} & (17) \\ {{{LNSR}_{M}\lbrack i\rbrack} = {\sum\limits_{j = 1}^{N}{{LNSR}_{j}\lbrack i\rbrack}}} & (18) \end{matrix}$

Equations (16) and (17) demonstrate the use of the chain rule in the partial derivative, and introduce the concept of a modem coefficient A[i] that encapsulates receiver modem information. The modem coefficient A[i] multiplies the partial derivative of the noise (LNSR) to power ratio of a specific section.

Through proper scaling of the metric D[i], the modem coefficient A[i] can be made equal to the first derivative of the example concave value function ƒ. A[i]=ƒ¹(SNR_(M)[i]). The receiver modem whose channel i traverses the section is capable of determining the value of the modem coefficient A[i]. For channels where receiver modem information is unavailable, the modem coefficient A[i] can be set to equal the number 1.

The value of this proper scaling of the metric D[i] is to ground the example concave value function ƒ of measured margin onto the example objective function V₁ given in Equation (1) which can be shown to either maximize capacity or the capacity-product depending on the choice of the weighting value C[i]. When the modem coefficient A[i] is set to equal the number 1 for all channels, the derivative in Equation (16) for the example objective function V₂ is identical to the derivative in Equation (4) for the example objective function V₁. Thus the example concave value function ƒ which tends to help channels with less margin at the expense of channels with more margin will operate around the control parameters that are close to either maximizing capacity or the capacity-distance product of the optical network.

The following discussion derives the gradients of the example objective function V₂, which are evaluated to provide a direction for adjustment of control parameters. A gradient vector ∇V_(2j) for control of the ingress WSS node 12 is derived. A gain gradient for TOP control of the optical amplifier devices 18 is derived.

By comparing Equation (16) and Equation (3), it is apparent that the gradients derived for the example objective function V₁ are applicable to the example objective function V₂, with the insertion of the modem coefficient A[i]. In cases where the modem coefficient AN equals 1 for all channels, the gradients derived for the example objective function V₁ are identical to the gradients derived for the example objective function V₂.

The partial derivative of the example function V₂ with respect to WSS loss for the channel kin the section is therefore given by Equation (19):

$\begin{matrix} {\frac{\partial V_{2}}{\partial{P\lbrack k\rbrack}} = {\sum\limits_{j = 1}^{N}\left\lbrack {{{- {A\lbrack i\rbrack}}{C\lbrack k\rbrack}\;\frac{{ASE}_{j}\lbrack k\rbrack}{e^{P_{j}{\lbrack k\rbrack}}}} + {\sum\limits_{k = 1}^{N_{CH}}{2\;{A\lbrack i\rbrack}{C\lbrack k\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}} & (19) \end{matrix}$ where A[i]=ƒ¹ (SNR_(M)[i]).

A gradient vector ∇V₂ for the section comprises the partial derivative ∂V₂/∂P[k] for each channel k from 1 to N_(CH). The gradient vector ∇V₂ can be evaluated from the modem coefficients A[i], the customer values C[i], the known characteristics NL_(j)[i,k], and the measured or estimated data {P_(j)[i]} and {ASE_(j)[i]}.

The gain gradient for the span j of optical fiber, averaged over all wavelengths, can be formed as given by Equation (20):

$\begin{matrix} {\frac{\partial V_{2j}}{\partial G_{j}} = {\sum\limits_{k = 1}^{N_{CH}}\left\lbrack {{{- {A\lbrack i\rbrack}}{C\lbrack k\rbrack}\frac{{ASE}_{j + 1}\lbrack k\rbrack}{e^{P_{j + 1}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{A\lbrack i\rbrack}{C\lbrack i\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}} & (20) \end{matrix}$ where A[i]=ƒ¹(SNR_(M)[i]).

The gain gradient for the optical amplifier device j, averaged over all wavelengths, can be evaluated from the modem coefficients A[i], the customer values C[i], the known characteristics NL_(j)[i,k], and the measured or estimated data {P_(j)[i]} and {ASE_(j)[i]}.

The values of the gradients are used in steepest descent algorithms to adjust control parameters of the section by a small step in a direction of optimization of the objective function.

Small adjustments are applied to loss values of a WSS node and to target TOP values of optical amplifier devices. The steepest descent algorithm is applied to the WSS node while assuming that the gains of the optical amplifier devices are fixed. The steepest descent algorithm is applied to all of the optical amplifier devices in parallel while assuming that the WSS pixel drive values are fixed.

For example, two loops may be run in parallel with a decoupling factor, as expressed in the vector Equation (21), Equation (22) and Equation (23):

$\begin{matrix} {{WSS\_ PowerTarget}_{NEW} = {{WSS\_ PowerTarget} - \left( {\frac{MAXSTEP}{\max\left( {\bigtriangledown\; V_{2}} \right)}*\left\lbrack {{\bigtriangledown\; V_{2}} - {{mean}\left( {\bigtriangledown\; V_{2}} \right)}} \right\rbrack} \right)}} & (21) \\ {{TOP\_ Target}_{NEW} = {{TOP\_ Target} - {{{sign}\left\lbrack \frac{\partial V_{2j}}{\partial G_{j}} \right\rbrack}*0.1*{MAXSTEP}}}} & (22) \\ {{{{if}\mspace{14mu}{TOP\_ Target}_{NEW}} \geq {TOP}_{LIMIT}},{{{set}\mspace{14mu}{TOP\_ Target}_{NEW}} = {TOP}_{LIMIT}}} & (23) \end{matrix}$ where WSS_PowerTarget_(NEW) and WSS_PowerTarget have values for each channel k from 1 to N_(CH), the decoupling factor in this example is 0.1, and the target TOP for the optical amplifier device j is subject to an upper limit. An example MAXSTEP is 0.2 dB.

TOP control is used to decouple incremental SNR optimization in this section from changes occurring other sections of the optical fiber network. The subtraction of the change in average power (which is denoted mean (∇V₂) in Equation (21) but is not quite equal to the average of the changes) keeps the WSS output power constant, and the power launched into the first span is controlled by the TOP of the first amplifier. This scaling is important when one or more of the amplifiers reaches the limit of their TOP and can provide no more power. With this scaling, the allocation of that limited power between the wavelengths is cleanly optimized.

Note also that there is no reliance on a second derivative for step size, and that this simple algorithm is robust to noise.

When the ability for rapid introduction of new channels is desired, idlers may be used to pre-allocate the effects of those channels.

An ASE idler is treated by the first aspect (example objective function V₁) as any other channel, with the appropriate XPM generator coefficient. The channel weight could be set to a static value of one. In the second aspect (example objective function V₂), once a modem signal is switched to replace this ASE, then the margin from that modem would be used to calculate the new weight in the usual way. The diminished default value is used again after the ASE is switched back in.

To not cause the XPM from ASE idlers, a limited number of virtual idlers can be used. Virtual idlers are treated just like ASE idlers, except that their XPM generator coefficient is set to equal that of the modulation expected to be used. Virtual idlers do not consume photons, so the TOP limits need to be reduced by the virtual wattage.

A Boolean acceptance criteria should be used to decide on the choice of a virtual idler versus an ASE idler in order to limit the SRS impact of their sudden conversion to real signals, assuming that all virtual idlers are allowed to switch at once. Define A to be the integral of virtual power spectral density out of the WSS, across a 1 THz interval centered at the wavelength of the candidate virtual idler, including the virtual power of that candidate idler. Define B to be the integral of real power spectral density across the same 1 THz interval centered at the wavelength of the candidate virtual idler. Choose epsilon to be a small positive number to avoid division by zero, e.g. 100 microWatts. The virtual idler is acceptable if A/(B+epsilon)<0.25.

In yet another aspect, the objective function is a combination of the above two objective functions. For example, the example objective function is given by Equation (24): V ₃ =V ₂ −V ₁=Σ_(i=1) ^(N) ^(CH) C[i]D[i]ƒ(SNR_(M) [i])−Σ_(i=1) ^(N) ^(CH) C[i]LNSR[i]  (24)

With this example objective function V₃, the goal of the objective function is to balance the goals of minimizing the total degradation through the section with maximizing capacity or reliability or both of the optical fiber network by re-allocating margin among the channels that are propagated through the section. The discussion above of deriving gradients and applying the gradients in steepest descent algorithms is applicable also to the example objective function V₃.

Returning now to FIG. 2, consider how this example section 10 could be modified to independently amplify different bands of transmission. For example, the section 10 could simultaneously handle the C band (1525 nm to 1565 nm) and the L band (1565 nm to 1610 nm). The ingress WSS node 12 could have two independent WSS filters to control individual channel powers for the C band and the L band, respectively. Each of the optical (pre-)amplifier devices 18,20 could be replaced by a set of two optical (pre-)amplifier devices, one for the C band and one for the L band. The C-band channels and the L-band channels propagate through the same spans 16 of optical fiber, where there is fiber nonlinear interaction between all the channels. That is, the nonlinear interaction in the spans of optical fiber is across all channels being propagated, including C-band channels and L-band channels. There is also strong SRS which makes for significant power differences between channels compared to the case of a single band, given that the SRS is approximately proportional to the square of the optical bandwidth.

In a variation of the first aspect, the example objective function V₁ applies to the full set of channels in the C band and the L band.

The partial derivative of the example objective function V₁ with respect to C-band WSS loss for the channel k in the section is given by Equation (25), where the channel k is in the C band:

$\begin{matrix} {\frac{\partial V_{1}}{\partial{P\lbrack k\rbrack}} = {{\sum\limits_{j = 1}^{N}\frac{\partial V_{1j}}{\partial{P_{j}\lbrack k\rbrack}}} = {\sum\limits_{j = 1}^{N}\left\lbrack {{{- {C\lbrack k\rbrack}}\;\frac{{ASE}_{j}\lbrack k\rbrack}{e^{P_{j}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{C\lbrack k\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}}} & (25) \end{matrix}$ where P[k] is the power of the channel k out of the ingress WSS node 12 which affects all spans in the section.

The partial derivative of the example objective function V₁ with respect to L-band WSS loss for the channel k in the section is given by Equation (26), where the channel k is in the L band:

$\begin{matrix} {\frac{\partial V_{1}}{\partial{P\lbrack k\rbrack}} = {{\sum\limits_{j = 1}^{N}\frac{\partial V_{1j}}{\partial{P_{j}\lbrack k\rbrack}}} = {\sum\limits_{j = 1}^{N}\left\lbrack {{{- {C\lbrack k\rbrack}}\;\frac{{ASE}_{j}\lbrack k\rbrack}{e^{P_{j}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{C\lbrack k\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}}} & (26) \end{matrix}$ where P[k] is the power of the channel k out of the ingress WSS node 12 which affects all spans in the section.

In Equation (25), the summation of the nonlinear interaction is over all N_(CH) channels i in the C band and in the L band. In Equation (26), the summation of the nonlinear interaction is over all N_(CH) channels i in the C band and in the L band.

A gradient vector ∇V₁(C) for the section for the C band comprises the partial derivative ∂V₁/∂P[k] for each channel k in the C band from 1 to N_(CH) ^(C). A gradient vector ∇V₁(L) for the section for the L band comprises the partial derivative ∂V₁/∂P[k] for each channel k in the L band from 1 to N_(CH) ^(L). The gradient vectors ∇V₁(C) and ∇V₁(L) can be evaluated from the customer values C[i], the known characteristics NL_(j)[i,k], and the measured or estimated data {P_(j)[i]} and {ASE_(j)[i]}.

The gain gradient for the optical amplifier device j, averaged over all wavelengths in the C band, can be formed as given by Equation (27):

$\begin{matrix} {\frac{\partial{V_{1j}(C)}}{\partial G_{j}} = {\sum\limits_{k = 1}^{N_{CH}^{C}}\left\lbrack {{{- {C\lbrack k\rbrack}}\frac{{ASE}_{j + 1}\lbrack k\rbrack}{e^{P_{j + 1}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{C\lbrack i\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}} & (27) \end{matrix}$ where the outer summation is over the channels k in the C band, and the inner summation is over all N_(CH) channels i in the C band and the L band.

The gain gradient for the optical amplifier device j, averaged over all wavelengths in the L band, can be formed as given by Equation (28):

$\begin{matrix} {\frac{\partial{V_{1j}(L)}}{\partial G_{j}} = {\sum\limits_{k = 1}^{N_{CH}^{L}}\left\lbrack {{{- {C\lbrack k\rbrack}}\frac{{ASE}_{j + 1}\lbrack k\rbrack}{e^{P_{j + 1}{\lbrack k\rbrack}}}} + {\sum\limits_{i = 1}^{N_{CH}}{2\;{C\lbrack i\rbrack}{{NL}_{j}\left\lbrack {i,k} \right\rbrack}e^{2\;{P_{j}{\lbrack i\rbrack}}}}}} \right\rbrack}} & (28) \end{matrix}$ where the outer summation is over the channels k in the L band, and the inner summation is over all N_(CH) channels i in the C band and the L band. The total number of channels in the C band and the L band, denoted N_(CH), is the sum of the number of channels in the C band, denoted N_(CH) ^(L), and the number of channels in the L band, denoted N_(CH) ^(L). That is, N_(CH)=N_(CH) ^(C)−N_(CH) ^(L).

The gain gradients can be evaluated from the customer values C[i], the known characteristics NL_(j)[i,k], and the measured or estimated data {P_(j)[i]} and {ASE_(j)[i]}.

Small adjustments are applied to loss values of a WSS node and to target TOP values of optical amplifier devices. The steepest descent algorithm is applied to the WSS node while assuming that the gains of the optical amplifier devices are fixed. The steepest descent algorithm is applied to all of the optical amplifier devices in parallel while assuming that the WSS pixel drive values are fixed.

For example, four loops may be run in parallel with a decoupling factor, as expressed in the vector Equations (29) and (30), Equations (31) and (32) and Equations (33) and (34):

$\begin{matrix} {{{WSS\_ PowerTarget}_{NEW}(C)} = {{{WSS\_ PowerTarget}(C)} - \left( {\frac{MAXSTEP}{\max\left( {\bigtriangledown\;{V_{1}(C)}} \right)}*\left\lbrack {{\bigtriangledown\;{V_{1}(C)}} - {{mean}\left( {\bigtriangledown\;{V_{1}(C)}} \right)}} \right\rbrack} \right)}} & (29) \\ {{{WSS\_ PowerTarget}_{NEW}(L)} = {{{WSS\_ PowerTarget}(L)} - \left( {\frac{MAXSTEP}{\max\left( {\bigtriangledown\;{V_{1}(L)}} \right)}*\left\lbrack {{\bigtriangledown\;{V_{1}(L)}} - {{mean}\left( {\bigtriangledown\;{V_{1}(L)}} \right)}} \right\rbrack} \right)}} & (30) \\ {{{TOP\_ Target}_{NEW}(C)} = {{{TOP\_ Target}(C)} - {{{sign}\left\lbrack \frac{\partial{V_{1j}(C)}}{\partial G_{j}} \right\rbrack}*0.1*{MAXSTEP}}}} & (31) \\ \begin{matrix} {\mspace{79mu}{{{{if}\mspace{14mu}{TOP\_ Target}_{NEW}(C)} \geq {{TOP}_{LIMIT}(C)}},}} \\ {\mspace{79mu}{{{set}\mspace{14mu}{TOP\_ Target}_{NEW}(C)} = {{TOP}_{LIMIT}(C)}}} \end{matrix} & (32) \\ {{{TOP\_ Target}_{NEW}(L)} = {{{TOP\_ Target}(L)} - {{{sign}\left\lbrack \frac{\partial{V_{1j}(L)}}{\partial G_{j}} \right\rbrack}*0.1*{MAXSTEP}}}} & (33) \\ \begin{matrix} {\mspace{79mu}{{{{if}\mspace{14mu}{TOP\_ Target}_{NEW}(L)} \geq {{TOP}_{LIMIT}(L)}},}} \\ {\mspace{79mu}{{{set}\mspace{14mu}{TOP\_ Target}_{NEW}(L)} = {{TOP}_{LIMIT}(L)}}} \end{matrix} & (34) \end{matrix}$ where WSS_PowerTarget_(NEW)(C) and WSS_PowerTarget(C) have values for each channel k in the C band from 1 to N_(CH) ^(C), WSS_PowerTarget_(NEW)(L) and WSS_PowerTarget(L) have values for each channel k in the L band from 1 to N_(CH) ^(L), the decoupling factor in this example is 0.1, and the target TOP for the optical amplifier device is subject to an upper limit (dependent on the band). An example MAXSTEP is 0.2 dB.

In a variation of the second aspect, the example objective function V₂ applies to the full set of channels in the C band and the L band. Similar equations and loops can be derived for that case, for example, by replacing the customer values C[i] in Equations (25) through (34) with the product of the customer values C[i] and the modem coefficients A[i].

For clarity, the examples apply a Gaussian nonlinearity noise model. The methods described in the document can be used where other models of optical nonlinear interactions provide a better representation, such as where the nonlinearities are not substantially independent between spans.

For clarity, the examples use ASE from Erbium doped fiber amplifiers (EDFAs). Other power or gain dependent degradations such as double-Raleigh scattering from Raman amplifiers may be included in the estimations.

Modern high capacity optical transmission systems use coherent modems (also known as coherent transceivers). The techniques described in this document may also be used with other kinds of optical transmitters and receivers.

The description shows specific examples of objective functions and the derivation of control algorithms from those objective functions. Other objective functions may be used. Using a convex objective function and deriving a control algorithm from an objective function is convenient. However, other control algorithms and methods may be used.

For simplicity, the control algorithms described in this document have sole control of the control parameters of the section. Other control algorithms or provisioning or constraints may also be present or active.

The scope of the claims should not be limited by the details set forth in the examples, but should be given the broadest interpretation consistent with the description as a whole. 

What is claimed is:
 1. A method for adjustment of one or more control parameters of a section of an optical fiber network, the method comprising: taking measurements of optical signals in the section; deriving estimated data from the measurements and from knowledge of the section, where the estimated data is a function of optical nonlinearity and of amplified spontaneous emission; evaluating gradients of an objective function using the measurements and the estimated data; and applying one or more control algorithms using at least the gradients to adjust the one or more control parameters.
 2. The method as recited in claim 1, wherein the control parameters include per-channel power control.
 3. The method as recited in claim 2, wherein the control parameters include loss values of a wavelength selective switch (WSS) node in the section.
 4. The method as recited in claim 1, wherein the control parameters include total power control.
 5. The method as recited in claim 4, wherein the control parameters include target total output power values or target gain values of one or more optical amplifier devices in the section.
 6. The method as recited in claim 1, wherein the control parameters include per-channel power control and total power control, and the control algorithm to adjust per-channel power control is substantially decoupled from the control algorithm to adjust total power control.
 7. The method as recited in claim 1, wherein the estimated data includes estimated per-channel optical power or estimated power spectral density at the output of one or more optical amplifier devices in the section.
 8. The method as recited in claim 1, wherein applying the one or more control algorithms to adjust the one or more control parameters comprises adjusting the one or more control parameters in a direction of one or more dimensions of those gradients.
 9. The method as recited in claim 8, wherein the optical signals comprise two or more bands of channels that are independently amplified through the section, and optimization of the objective function jointly optimizes over all channels of the two or more bands.
 10. The method as recited in claim 8, wherein a goal of the objective function is to minimize degradations through the section.
 11. The method as recited in claim 10, wherein optimization of the objective function minimizes a weighted sum of ratios of the total noise power from amplified spontaneous emission and optical nonlinearities to the power of the optical signals.
 12. The method as recited in claim 11, wherein optimization of the objective function maximizes a capacity-distance product of the optical fiber network.
 13. The method as recited in claim 11, wherein optimization of the objective function maximizes capacity of the optical fiber network.
 14. The method as recited in claim 8, wherein a goal of the objective function is to maximize capacity or reliability or both of the optical fiber network by reallocating margin between channels.
 15. The method as recited in claim 14, wherein the objective function incorporates a concave value function of margin of one or more channels propagated through the section.
 16. The method as recited in claim 15, wherein for a channel propagated through the section, the gradients incorporate a function of information provided by a receiver modem that receives that channel.
 17. A method for adjustment of one or more control parameters of a section of an optical fiber network, the method comprising: taking measurements of optical signals in the section; deriving estimated data from the measurements and from knowledge of the section, where the estimated data is a function of optical nonlinearity and of amplified spontaneous emission; forming an objective function that combines performance of multiple channels; and applying one or more control algorithms using the estimated data and the objective function to adjust the one or more control parameters.
 18. The method as recited in claim 17, wherein the measurements include a total output power.
 19. The method as recited in claim 17, wherein the measurements include per-channel optical power or power spectral density.
 20. The method as recited in claim 19, wherein the measurements of per-channel optical power or power spectral density are taken at fewer geographic locations than a number of geographic sites in the section having optical amplifier devices.
 21. The method as recited in claim 20, wherein measurements of per-channel optical power or power spectral density are taken at two geographic locations and the section has more than three geographic sites having optical amplifier devices.
 22. The method as recited in claim 17, wherein the optical fiber network comprises multiple sections, and multiple instances of the one or more control algorithms operate in parallel.
 23. The method as recited in claim 22, wherein a topology of the optical fiber network comprises a branch topology.
 24. The method as recited in claim 22, wherein a topology of the optical fiber network comprises a mesh topology. 